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Abstract —We address the problem of uplink co-operative 
reception with constraints on both backhaul bandwidth and the 
receiver aperture, or number of antenna signals that can be 
processed. The problem is cast as a network utility (weighted 
sum rate) maximization subject to computational complexity and 
architectural bandwidth sharing constraints. We show that a 
relaxed version of the problem is convex, and can be solved 
via a dual-decomposition. The proposed solution is distributed 
in that each cell broadcasts a set of demand prices based on 
the data sharing requests they receive. Given the demand prices, 
the algorithm determines an antenna/cell ordering and antenna- 
selection for each scheduled user in a cell. This algorithm, 
referred to as LiquidMAAS, iterates between the preceding two 
steps. Simulations of realistic network scenarios show that the 
algorithm exhibits fast convergence even for systems with large 
number of cells. 

I. Introduction 

Uplink Co-ordinated Multi-Point (CoMP) is a promising 
technique for increasing the capacity of 4G networks (T), 
0, 0, SI- The uplink is gaining increasing attention due 
to the dramatic increase of user-generated data in the form of 
photos, videos and file-sharing. In practice, sharing of uplink 
received signals across cells is limited by backhaul bandwidth. 
In addition, the receiver aperture, or number of signals from 
antennas at neighboring cell sites that can be processed at a 
particular cell, may be limited due to hardware constraints. 

Several approaches for signal sharing and combining have 
been proposed for uplink CoMP. A performance analysis of 
different combining methods with different backhaul band¬ 
width requirements is presented in 0. To reduce the amount 
of sharing, dynamic clustering of cells has been proposed, e.g., 
in 0, 0- In addition to limiting the amount of information 
that can be shared across cells, the limited backhaul bandwidth 
also introduces latency, which is addressed in |8l . 

Previous work on CoMP has generally assumed that the set 
of cells that share information is fixed a priori by the topology 
and the bandwidth of the backhaul links. In this paper we relax 
this assumption and consider the problem of optimizing sets of 
helper cells that pass along their uplink signals to other cells. 
(A cell can both share its signals as a helper cell and receive 
signals, or help from other cells.) We account for architectural 
constraints that limit the set of potential helper cells (which 
differs from cell to cell), backhaul constraints that limit the 
number of cells to which a helper cell can send its received 
signals {egress constraint), and hardware constraints at the cell 


site, which may limit the number of incoming signals the cell 
can combine (ingress constraint). 

Related work in 0, 0 has considered the problem of 
MMSE receiver estimation under compression and backhaul 
constraints. We focus here instead on a simpler sharing for¬ 
mulation that introduces explicit constraints on egress band¬ 
width, to arrive at a convex weighted sum rate formulation 
with guaranteed convergence. Imposing these explicit egress 
bandwidth and ingress aperture constraints captures impor¬ 
tant architectural limitations in centralized or distributed co¬ 
operative networks. Our problem is then to maximize a sum 
rate objective subject to these architectural constraints. Since 
the cell sites are assumed to have multiple antennas, we refer 
to this problem as Multi-Antenna Aperture Selection (MAAS) 
for joint reception (IR)-CoMP (see also flOl ). 

The optimization of sets of helper cells is an integer 
program. Assuming max-ratio combining of received signals 
and relaxing the integer constraints, the optimization problem 
becomes convex. We present a distributed algorithm in which 
each helper cell announces an egress price, indicating the 
demand for its signals to help other cells, and each assisted cell 
computes an ingress price, indicating the potential improve¬ 
ment from adding a helper cell. The prices are used to compute 
an ordering of users/cells, which is then used to allocate 
the helper antennas across assisted cells. This assignment is 
iterated with updates for the prices. 

We refer to this algorithm as LiquidMAAS due to its ability 
to flexibly allocate help based on network load conditions, and 
show that it converges to the optimal allocation of helper cells 
across the network. Numerical results are presented that show 
that the algorithm converges in relatively few iterations even 
for a system with a large number of cells. Furthermore, the 
gains relative to an a priori fixed allocation of helper cells can 
be substantial. 

II. Multi-Antenna Aperture Selection 
A. Problem Setup 

Consider a network of J cells, each with multiple antennas. 
Each cell j serves a set of users u(j). The uplink signal from 
each user k £ u(j) is typically strongest at its serving cell; 
however, the signal could also be received with significant 
strength at other cells depending on the user location and 
network topology. We assume the antenna-combined signal 
for a particular user k £ u(j) can be shared with other cells. 


and is, of course, provided to that user’s serving cell for CoMP 
combining. 

Figure CD shows an example of a network with uplink 
data sharing. There are four cells with seven users, where 
u(l) = {1, 2}, u( 2) = {3}, u( 3) = {4, 5, 6}, and u( 4) = {7}. 
Let Afii(k) denote the set of cells with data that can be 
requested by the cell serving user k (ingress neighborhood). 
That typically corresponds to the set of cells where the user’s 
SINR is above a minimum threshold value (usually -10 dB). 
That is, J\f R (k) = {i G J,i ^ a(k) : > S'min}, where 

S^ 3 is the SINR at cell i for a user k in cell j (i.e., k G 
and er(fc) is the serving cell for user k. Note that A f R {k) may 
not be the actual set of helper cells for user k due to backhaul 
and aperture constraints. 

Also shown in Figure ([[} are the possible sharing variables 
along with the inter-connect. Note that cells 1 and 3, and cells 
2 and 4 do not exchange information. The sharing variable 
G {0,1} indicates whether or not cell i shares data 
for user k in cell j over the backhaul link * — > j. For the 
example shown in Figure (jTJ. the local neighborhood sets for 
the different users are 



Backhaul linkGall->Gall 


Fig. 1. Example cell topology with backhaul inter-connects. 


B. Network Utility Maximization 



Me(l)=Mi(2) = {2,4}; Affl(3) = {1,3}; 
■A/it(4) = A/}t(5) = A/it(6) = {2,4}; Mr(7) = {1, 3} 

The maximum receiver aperture size, or number of helper 
cells, for user k is denoted as L R {k), and is the number 
of cells in the neighborhood JV R (k). For a cell i, its egress 
neighborhood to cell j is denoted by A'r(*, j), and is the set 
of users in cell j it can potentially help, ignoring backhaul 
and aperture constraints. That is: 

A/r(f, j) = {k G u(j) : i G U R {k)} ; i,j G J, i ± j 

For the example in Figure ([[}. L R {k) = 2 V/c G u(j), Vj = 
{1, 2, 3,4}. The egress neighborhoods are: 

A/t(1j 2) = {3}; A/t(1,4) = {7}; A/'t (2, 1) = {1,2}; 
A/t( 2, 3) = {4,, 5,6}; A/r(3,2) = {3}; A/'t (3,4) = {4,5,6}; 

At(4,1) = {1,2}; A/t(4, 3) = {4, 5,6}; 

Note that if there is no possibility of sharing between two 
cells, then the corresponding egress set is null, i.e., A/r(2,4) = 
A/t( 4, 2) = 0. The notation is summarized in Table [Q 

TABLE I 

Variable Definitions 


Our problem is to select the set of helper cells for each 
user k to maximize an overall network objective, subject to the 
egress (backhaul bandwidth) and ingress (aperture) constraints. 
This selection is assumed to occur for a fixed schedule of users 
in a given transmission time, and after reception of the signals. 
We assume that a user is always processed by its serving 
cell’s antennas. If another cell’s antennas are included in the 
aperture (helper set) for a user, then that cell performs local 
processing (e.g., Max-Ratio Combining (MRC) or Minimum 
Mean Squared Error combining across its own antennas), and 
forwards those post-combined signals to the serving cell. That 
cell in turn performs MRC of the signals from different cells. 
Therefore, the combined SINR is the sum of the SINRs from 
the serving and helper cells. 

We therefore obtain the following optimization problem 
(NUM) over the sharing variables 


Maximize 


E E UJkPk log 

j&J kGu(j) 


1 + sjUi 


E 

*€.A/jl(fc) 


Qk k 


Subject to: 


(i) 


J 

Set of all cells in the network, with J = \J"\. 

u ti) 

Set of users connected to cell j. 

cr(k) 

Serving cell for user k. 

X i^j 

BW fraction that cell i shares with user k in cell j. 

cAc 

SINR experienced by user k at cell i. 

Pk 

Bandwidth fraction allocated to user k. 

uk 

Scheduling priority weight for user k. 

La 

Maximum number of helper cells allowed. 

Afl(fc) 

Ingress neighborhood for user k in cell j. 


Maximum aperture of helper cells (= |Af^(/c)|). 

At(l j) 

User egress neighborhood for cell i to cell j. 

l/rjy 

Maximum per-cell egress bandwidth. 


E x i-rj - La V -7 G Vfc G u (j ) ( 2 ) 

ieN R (k) 

Y, E A xL] < Lf Vi G J (3) 

k£Ar T (i,j) 

e 1 J’ V -7 G 31 Vfc e u U)>^ j 14) 

where wt and ffr are defined in Table U 

The cost function is concave and the inequality constraints 
are linear. If we relax the variables to be 0 < < 1, 

the optimization problem is convex. A fractional value for the 
sharing variable has an interpretation of a helper cell sending 























only part of the resource block (RB) allocated to a user. Since 
this is acceptable within the inter-connect architecture, and 
also for the receiver processing, we relax the variables to 
obtain a convex formulation. In fact, we will show that in the 
solution to the relaxed problem, for each user k, all sharing 
variables except one satisfy the integer constraint. 

The ingress constraint is a computational complexity con¬ 
straint based on the maximum number of cell signals that can 
be processed by the receiver, and is independent of the allo¬ 
cated bandwidth fraction of the useiQ. The egress constraint, 
on the other hand, is a backhaul bandwidth constraint, and 
hence accounts for the fraction of bandwidth a user’s signal 
occupies. 

Note that the egress constraint © is the only coupling 
constraint. If that were absent, the solution would be straight¬ 
forward: for each user k £ u(j), order the cells according to 
the SINRs for that user, and pick the min{L^ 4 , L R {k)} top 
cells for receiver processing. 

III. Characterization of Solution 

We now present a distributed algorithm for solving the 
preceding NUM problem. Since the relaxed problem is convex, 
we seek a set of sharing variables that satisfies the KKT 
conditions. Define a combined SNR metric for user k as 


The KKT conditions for the NUM are then given by: 


uikPkS^j 

g(*k) 


X k 'fiftk 7a T @i,k 



£ £ 

l$Lj raEA/r(M) 


= o 

= o 


( 8 ) 

(9) 


A k\L A - £ x n-tj = 0 (10) 

\ n£Nn(k) ) 

(i — x i^j) = o (11) 

8i,kxLj = 0 ( 12 ) 


for every (i. j) £ J,i j,k £ u(j). Thus the NUM 
decomposes into a set of local per-user optimization problems 
for a given a set of egress prices {ipi}. 

We can now apply a dual decomposition to obtain a dis¬ 
tributed iterative algorithm for solving the NUM. Namely, 
given set of shadow prices {ipi(t)} and {Afc(£)} for the egress 
and ingress constraints (1312b at iteration t, each user solves 
its cell selection problem locally by maximizing the objective 
C k (x k , A k(t), 4'■(£), Tfc, 0fc). Define this primal solution to be 
{x£(£), T£(£), ©£(£)}. Then, the egress and ingress prices can 
be updated using the sub-gradient algorithm 


< 7 (x fc ) = 1 + SjLtf + £ S k m ^x k m ^ (5) 

m£A/B(fc) 

where x*. = {x k ^j} i( zjj- R ( k y We can then wf ite the La- 
grangian for the NUM as 

£(x,A,^,r,e) = £ £ c fc (x fc , A fc ,*,r fc ,e fc ) 

j£J keu(j) 

+ £ ipi Lf (6) 

iej 

where the cost function for user k is 

£fc(xfc, A*,, T, T k , Ok) = Wfc/3fc log [ff(x fc )] 

+Afe I L a - £ J -0 k £ ipi 
\ f R (k) J i£Nn(k) 

+ £ 7 a (i-*££+ £ V) 

i&J^Rik) i£NR.(k) 

the Lagrange multipliers ipi, X k ,'Yi.fc, 0i,k > 0, {T k , 0 fe } is the 
set of KKT multipliers per user, and {'T, A, T, 0} is the entire 
set of KKT multipliers. Note that to write the Lagrangian C 
this way, the egress constraint must be broken into partial sums 
of sharing variables for a particular user k: 


ipi{t + 1) = [ipi{t) - vAi(t)\_ 


A i(t) = Uf- £ £ ^C(i) 13) 

y /£ m€.Afr(i,l) / 

X k (t+l) = [X k (t)-v(L R (k)- £ x k n % 3 m+ (14) 

n^Nn(k) 


where v > 0 , is a small step-size, and [y] + = max{y, 0 }. 

It remains to solve the primal problem in ©. For this define 
the metric: 

Pk + A k ,, 

= - a ok ■ d 5 ) 

UJ k Pk 


Lemma 3.1: If ^ rni k , VZ ^ i, uj k > 0, 0 k > 0, 

and .a > 0 Vfc, then there can be at most one variable 
in that has a fractional value. All others are 

either 0 or 1 . 

Proof: Consider the ingress problem for user k £ u(j), 
that is, maximization of the objective ©. Assuming a fixed 
set of egress prices {ipi}, i £ Af R {k ), we rewrite the KKT 
condition ® as 


ff(xfc) 


— tTU k T 


V* £ Af R {k) since u} k P k S k ^j > 0, and where 


(16) 


£^ £ £ faxi^j = £ £ p* £ 

iej ZcEA/t(«J) jej keu(j) i£Afii(k) 


$i,k 


@i,k 

uikPkSi^j 


( 17 ) 


*Even ignoring computational limitations, limiting the number of signals 
that are combined can improve performance due to finite training, associated 
channel estimation error, and receiver imperfections na- 


it follows from the KKT conditions that 

• if a£j- = 1 , then 7 i)fc > 0 , 0 i)fc = 0 , and d i>k > 0 ; 

• if = 0, then 7 i]fc = 0, 6^ k > 0, and S t , k < 0; 






• if G (0,1), then j i>k = d itk = $ i)fe = 0. 

There cannot be two sharing variables and l 

that both take fractional values in the optimal solution, since 
(DU cannot be satisfied for both i and l with 5i, k = 6 i k = 0 
given that m iik ± mi, k V l / i; i,l G Af R {k). ■ 

Next we characterize the solution to the primal problem for 
user k. 

Definition Let x k be the optimal solution to i(7}. Define the 
set of active variables in this solution to be Z acl j ve , and the 
index of the cell with the maximum metric in this set as i a : 

^active = {t€Mt(fc):0<®^.<l} 

i Q = arg max nrii, k (18) 

ieI active 

Theorem 3.2: For a fixed X k and ipi, i G J\f R (k), assume 
w.l.o.g. that rrii^. {i = 1, 2,3,..., L R (k)}, is increasing for 
user k. Then, the sharing variables that maximize C k in {7} 
are given by: 

• x^j = 1 for i <i a , 

• Xil+j = 0 for i > i a , 

• G (0,1) if 1 +gfe^ — gn —- > m iotk , else 

r k * . — I 

Proof: Let the optimal shadow price difference be { 6 * k } 
corresponding to the optimal solution x* k . We show that 
-7-active {1; 2, , i 0 }■ 

PART P. 2i ac tive — {1 > 2, • ■ ■, *o}- 

Since mi, k is increasing in {* = 1,2,3 ,...,L R (k)}, it 
follows that to satisfy (IT6l) . 

6 *i, k > 6 * 2 , k >■■■> 5* ok >...> 6 * LR{k)>k (19) 

Since i a G I act j ve , 5f k > 0. Hence, from (IThl i. 

6* k >0 => x*^ = l Vi = l,2,...,* 0 -1 

=> {1,2, C I active (20) 

PART IP. I active C {1,2 

For this, we show that {i a + l, * 0 + 2,..., L R (k)} C I active , 
that is, x^j = 0 Vi = i 0 + 1,... L R {k). This follows by 
contradiction. Assume that i' G I act j ve , and i' > %„. By 
definition then 61, k > 0. Since i' > i a , we have S*, k < 6* o k 
from (fl9l) . This implies rrii’, k > m.i 0 , k , which contradicts the 
definition of i a in ©. 

From these, it follows that: 

-7-active = {f > 2> • • • > *o} 


Further, the one variable that can have 


value 


i+sjLi+IA 


r *k 

1 


and 


sE~+sf 


found as a 

T Te = Wlio.k’ 


a fractional 
solution to 


Remark The value of i 0 can be then found by a search such 
that the KKT criterion in (fl6l > is met with 5* k > 0 for i < i a , 
6* k < 0 for i > i 0 and S* o k > 0 with some 0 < < 1. 


IV. LiquidMAAS Algorithm 

The following algorithm is based on the preceding proper¬ 
ties of the solution. 

1) Initialize egress (bandwidth demand) and ingress 
(aperture) prices: 

a) tpi = e > 0 \/i 

b) A/c = 0 VI" G u(j), j G J. 

2) For each k G u(j), j G J , do: 

a) Order helper cells: For user k, order the helper 
cells in increasing value of metric rrii, k . Let that 
ordered list be l k = {ii,i 2 , ■ ■ ■,iL R (k )}• 

b) Select aperture: Initialize a;J • = 0 for i = 
l,--- ,L R (k). Add helper cells sequentially, i.e., 
for n = 1,2, • ■ • add a new helper cell n by setting 
x„_>j = 1, and check if 


1 


1 + 5: 






ok 

°i p ^j 


> m in , k (21) 


If true, then increment the aperture size to n. If 
false, then stop incrementing n and set 


where 


V = 


™i n ,k 


xt-H = max {0, y} 


- [1 + + J2l<n -1 ^ 


( 22 ) 




ok 


and y < 1, corresponding to the fractional sharing 
variable. 

3) Update Ingress Price: 

A fc t— max (Afc - v[L R {k) - L R (k)],0j 

where L R (k) = J2ieAr R (k) and ^ > 0 is a small 
step-size. 

4) Update demand headroom: Each cell computes the 
egress demand headroom: 


Aj — Lj’ 


E E 


(23) 


5) Update demand prices: 

V’i G- ma x(ipi - v\i, 0), 


6) Iterate: Repeat steps (2)-(5) until convergence (i.e., 
all | L R (k) — L R (k) \ < E\ and |Aj| < e 2 for some 
sufficiently small £i,e 2 > 0). 

Given the tolerances Ei,E 2 —> 0, if rrii, k ^ mi, k \/l f k , 
Wfc > 0, d k > 0, 5jE > 0 Vfc, then the LiquidMAAS algo¬ 
rithm is guaranteed to converge, and the allocation approaches 
the solution to the NUM problem in GB as £ 1,£2 0. 

This follows from the results in Lemma 13.11 Theorem 13.21 
and the subsequent remarks. Furthermore, we note that the 
algorithm can identify the integer-valued sharing variables in 
the solution within a finite number of steps (provided that 
£i, £2 are sufficiently small). The algorithm can be modified 









so that this optimality property also applies to the case where 
Wz.fc = for some l ^ k; however, we omit the details 
here. 

The algorithm is referred to as LiquidMAAS since it dis¬ 
tributes the helper cell load across the network, ensuring that 
all cells are able to meet their egress bandwidth constraint, 
while maximizing network utility within the ingress (aperture) 
constraint. In contrast, if we do the traditional cell selection 
approach that considers only SINRs at different cells, there 
could be situations where the demand on certain cells could 
be very high, much exceeding the egress bandwidth. If a cell 
is in high demand i.e., many cells request its antenna data, 
that cell’s ipi value will be high, indicating that there is a high 
price to be paid for getting that cell’s antenna data. If there is 
little demand for a cell i’s data, then ipi ~ 0. 

If all cells have the same demand, then the cell selection 
amounts to ordering purely by SINR. If the ingress price A/,. 
ipi Vi, then the ordering is approximately according to SINR. 
This ordering takes the potential help one can get from a helper 
cell in terms of SINR (S^. •) and weighs it against the price 
of obtaining it (ipi) and how tight the ingress constraint is 
(Afc). So if we have two cells that give about the same help 
in terms of SINR, one would pick the cell that is less “busy” 
or in demand (smaller ipi). If a cell is heavily in “demand”, 
some cells may choose to request from other cells that are in 
less demand, while still getting a reasonable improvement in 
SINR. 

The algorithm also automatically takes a user’s priority 
(wfc) into consideration while making a decision on the user’s 
aperture. This priority usually arises out of enforcing fairness 
for the different users via a utility such as proportional fairness. 
If a user is deemed to be high priority (large utk), then the 
aperture selection in (PH becomes small for that user, biasing 
the selection of larger apertures for such users. This is in 
agreement with what we expect intuitively, i.e., that high 
priority users should be able to pick larger apertures so as 
to obtain improved throughput. 

V. Simulation Results 

We evaluated the LiquidMAAS algorithm for the case of a 
standard 57-cell layout with an average of 10 users per cell. 
The simulation parameters are shown in Table fTil 

Users were allocated bandwidth equally based on the 
number of users connected to a given cell. A user’s SINR 
to all cells in the network was calculated based on this 
bandwidth allocation, and open-loop fractional power control 
EQ was assumed. The update for ipi and A& used a step- 
size of v = 0.005. Four algorithms were evaluated: (1) No 
CoMP, (2) MAAS without egress constraints, (3) MAAS with 
a randomized egress bandwidth control mechanism, and (4) 
LiquidMAAS. In case (2), each cell requests help based on 
their aperture limit, and helper cells grant any requested help. 
In (3), helper cells randomly grant requests for help until their 
egress constraint limit is reached. 

The ingress aperture constraint limit was set to La = 3, so 
that each user being gets a maximum of 3 helper cells’ data. 


TABLE II 

List of Simulation Parameters 


Simulation Parameter 

Values 

UL system bandwidth 

10 MHz 

Carrier frequency 

2.0 GHz 

Load average 

10 UE per cell 

Inter-site distance 

100m 

Receiver 

MRC 

Number of antennas per cell 

2 

Power control (a, Po ) 

(0.8, -80 dBm) 

Max transmit power 

24 dBm 

Noise figure 

4 dB 

Traffic model 

Full buffer 


The egress bandwidth limit L T was varied to see the effect 
on performance and convergence. The top part of Figure (0 
shows the convergence behavior of all cells’ egress bandwidth 
demand (for L T = 1.0), from which we observe that the 
algorithm converges in « 50 iterations. The convergence time 
is dependent on the step-size for updating the prices. It was 
observed that if the aperture and egress limits are similar, the 
step-size should be reduced to enable smoother convergence. 
The bottom part of Figure © shows the distribution of 
weighted sum rate (WSR) gain obtained by the various CoMP 
approaches. This shows that LiquidMAAS gives substantial 
gain over ’No CoMP’ and also over the randomized egress 
bandwidth control strategy. 

Figure © shows the distribution of the egress bandwidth 
and ingress apertures at convergence. We observe that the 
LiquidMAAS algorithm maintains tight control over the egress 
bandwidth, and adapts the ingress aperture to ensure that the 
constraints are met. Finally, Figure © shows the average WSR 
gain (over ’No CoMP’) as a function of the egress limit ( L T ). 
LiquidMAAS gives significant gains across the range of egress 
limits compared to the randomized strategy, and converges to 
the unconstrained MAAS solution for larger values of L T . 

VI. Conclusions 

In this paper, we presented an approach for joint algorithm 
and architecture optimization for co-operative communication 
networks. We considered the problem of uplink joint reception 
CoMP with backhaul bandwidth constraints, and showed that 
the resulting network utility maximization problem is convex. 
A distributed algorithm to solve this problem was presented, 
which consisted of local nodes computing their desired helper 
requests, followed by helper cells computing and updating an 
egress price for their bandwidth. At convergence, these itera¬ 
tions give a set of connections between helper and recipient 
cells, such that egress and ingress constraints are met at all 
cells, and overall network utility is maximized. Simulation 
results for a 57-cell system were presented to illustrate the 
efficacy of the approach, and show that the distributed dual 
decomposition-based gradient algorithm converges within tens 
of iterations. 








Demand convergence for all cells, Egress demand max =1 



Weighted Sum Rate Gain from Liquid MAAS and MAAS over No CoMF 



WSR Gain Per Cell (%) 


Fig. 2. MAAS algorithm performance comparison and LiquidMAAS egress 
bandwidth convergence for Lj, = 1.0. 
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Fig. 3. LiquidMAAS ingress aperture and egress bandwidth distribution for 
L t = 1.0. 
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Fig. 4. Algorithm performance comparison vs. egress bandwidth limit. 
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